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A SIMPLIFIED METHOD FOR ASSESSING THE EFFECT OF STEADY 

ROLLING ON ANGIE OF ATTACK AND SIDESLIP 

By Stanley F. Schmidt, Norman R. Bergrun, 

Robert B. Merrick, and Howard F. Matthews 


SUMMARY 


Presented herein is a simple method for analyzing the effects of 
inertia and aerodynamic cross -coupling on the response of airplanes in 
rolling maneuvers. The method is based on the concept that variations 
in the calculated steady-state angle of attack and sideslip with roll 
rate give a measure of the onset and degree of roll coupling. Compari- 
sons are made with the results from five-degree-of -freedom analog- computer 
studies of unaugmented and augmented airplanes , and it is shown that the 
method is useful in predicting trends but does not predict magnitudes with 
sufficient accuracy for loads analysis. Limited flight data are included 
also to demonstrate the use of the steady-state method in indicating 
unsatisfactory regions of flight. In addition, the method is shown to 
be useful in predicting the magnitude of the aileron deflection where 
serious roll coupling will occur. 


INTRODUCTION 


Recently, a number of high-speed airplanes have experienced large 
excursions in angle of attack and sideslip during a rolling maneuver 
(refs. 1 and 2) which have been attributed to inertia coupling. The 
possibility of such an occurrence was predicted in I 9 W by Phillips in 
reference 3 , wherein he showed that in a steady roll a divergence in yaw 
and pitch will occur when the rolling frequency approximately equals the 
uncoupled yawing or pitching natural frequencies . Although Phillips 
correctly defined the fundamentals of the inertia-coupling phenomenon, 
the development of a more comprehensive but still simple method for the 
prediction of serious roll-coupling effects was considered desirable. 

The purpose of this paper is to present such a method, referred to 
as the "steady- state method." This method is based on the concept that 
the variation of the steady-state angle of attack, sideslip, and other 
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quantities of* interest for various steady rolling velocities gives a 
measure of the onset and degree of roll coupling. The main substance of 
this paper , then, includes the derivation of formulas for the steady- 
state angle of attack, sideslip, etc., from the equations of motion, and 
an evaluation of this concept by means of a comparison of the steady- 
state results with analog- computer and flight data. 

Prior to the completion of this study, it had come to the authors 1 
attention that others (refs. 4, 5 , and 6 ) had also considered the steady- 
state approach. However, it is believed that this study treats the prob- 
lem more extensively and from a somewhat different point of view and 
therefore should be a useful addition to the expanding literature on the 
subject. 


NOTATION 


The analysis made in this report employs a system of body axes. The 
axis system is shown in figure 1. All angles are in radians unless other- 
wise noted. A dot above a symbol represents the first derivative with 
respect to time. 
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cyclic frequency in pitch and yaw, respectively, during rolling 
motion, cps 

acceleration due to gravity, ft/sec 2 
elevator and rudder servo gains, radians /volt 
altitude, ft 
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i x ,i Y ^z moment of inertia in roll, pitch, and yaw, respectively, 
^ referred to body axis, slug-ft 2 
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product of inertia referred to body axes, slug-f 


polar moment of inertia of jet engine rotating member, slug-ft* 
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mass of airplane, slugs 
Mach number 

rolling, pitching, and yawing velocities, respectively, 
radians /sec 

steady rolling velocity, radians/sec 
roll subsidence root, dimensionless 
wing area, sq ft 
time , sec 

incremental time interval, sec 
free-stream velocity, ft/sec 

voltage proportional to subscript quantity 

angle of attack 
angle of sideslip 

total aileron deflection, positive for trailing edge of right 
aileron down 

horizontal stabilizer deflection, positive for trailing edge 
down 

horizontal stabilizer deflection above that required for trim 
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rudder deflection, positive for left rudder 

declination of principal axis below body axis 

damping ratio in pitch and yaw, respectively, during rolling 
motion unless otherwise noted, dimensionless 

angle of pitch of longitudinal body axis 

air density, slugs /cu ft 

angle of roll 
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incremental angle of roll 

angular velocity of jet engine rotating member , radians/sec; 
positive rotation is clockwise vhen vieved from tail of 
aircraft 
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Lp^L^Iy^Lp &LJ2. (c lp , c Z 8a 8 a+Ci 5r &r, C l r 2V' % 2v) 

where Y,Z,L,M, and N are total forces or moments at the start of the 
rolling maneuver. 


THE STEADY- STATE METHOD 


Illustrated in figure 2 is a typical aircraft angle of attack versus 
sideslip analog- computer record for an applied constant aileron deflection 
with rudder and horizontal stabilizer held fixed at their initial trim 
values. In this instance, the aircraft started its rolling motion from 
initial trim conditions of 4° angle of attack and zero sideslip and when 
allowed to roll indefinitely tended to stabilize around a new steady-state 
condition, denoted in the figure by Og 8 and 0 8S . It is the variation of 
these and other pertinent steady-state values for various constant roll 
rates and the interpretation thereof that is the basis of the steady-state 
method which follows. - - 


Derivation of the Steady-State Formulas 


For this analysis a constant forward velocity and a constant altitude 
are assumed. The equations of motion used are the fiye-degree-of-freedcm. 


rigid-body equations . The five equations are: 

0+r-pa = Yp0+Y r r+YpP+Y g cos Q sin q>+Y (l) 

a-q+p0 = Z^a+ZgCOs Q cos cp+Z (2) 

q-I 1 pr+I 2 (p 2 -r 2 )+I I- fr = M^oM-M^oH-M^q+Mp 0 +M ( 3 ) 

r+IsPq+I^C-P+qrJ-Ijjq = N p 0+N^p+N r r+N p p+N (4) 

P+I 5 qr-I s (r+pq) = Lpp+LpP+Lpr+L (5) 


There are several different methods 1 for solving these equations to 
obtain the steady-state values for a constant roll rate (p = p Q ) and fixed 

■'■Another method is to transform equations (l) to (4) directly to 
algebraic equations by setting q = r= d= 0= Z=Yg = O and p = p Q . 
However, this procedure results in an equation which is of greater than 
the first degree in the particular quantity desired. For example, Oq S 
must be obtained from a cubic algebraic equation of which only one root 
is of practical significance. ^ _ 
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controls. A convenient one, which, has the advantage of yielding formulas 
(in the absence of nonlinear aerodynamic derivatives) that can be -used in 
some instances to easily predict effects of parameter cha n ges , is to 
linearize equations (l) to ( 4 ) by assuming: 


p = p 0 = a constant 

<p 


p dt = p Q t 


r 2 « p Q 2 


I^C[. ^1* 

0 a o° 


The resulting four linear differential equations then can be solved by 
standard operational techniques to give: 


(D 4 +a 3 D 3 +a 2 D 2 +a 1 D+a 0 )a = bo+b-jZgSin p Q t+b 2 Z g cos p Q t 

(D 4 +a 3 D 3 +a 2 D 2 +a 1 D+a 0 ) 3 = c 0 +c x Y g sin p 0 t+c 2 Y g cos p Q t 

where the a's, b’s, and c*s are constants and combinations of inertia 
and aerodynamic parameters and the rolling velocity, p Q . If these fourth- 
order differential equations are stable and the initial conditions are 
known, the time solution 2 for a and (3 will be of the form 


\ _ ^o (decaying (sinusoidal oscillations 

a ' ’ a D transients) due to gravity term Z g ) 

eft) = — + (decaying + (sinusoidal oscillations 
' a Q transients) due to gravity term Y g ) 


or the steady state is 


a. 


SS Sq 


Pss 


( 6 ) 

( 7 ) 


In reference 7 Sternfield has derived expressions of an approximate 
time solution for a and 3 ' assuming a constant roll rate. 
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The complete expressions for the constants a^, b Q , and c G may be found 
in Appendix A. The divergent boundary is given by a Q = 0, With appro- 
priate simplifications to the terms of ao, this equation can be reduced 
to that given by Phillips in reference 3 * 

To demonstrate certain, points in subsequent portions of this report, 
it is convenient to omit the minor parameters Y, Y p , Y p , Y r , Z, Np, and 
Mp in the expressions for a. Q , t> 0 , and c Q . The equations for Ogg and 
P gs then reduce to: 


(I 2 Pq 2 -M) (l 3 P 0 2 -%P 0 -Ng)-NM qPo 
(1 1 Po 2 -I M Po +I ^,)( I 3Po 2 - I #o' W 3 ) +M q N rPo 2 + M qZ a . N p 


( 8 ) 


(I 2 P 0 [ Cl3^a + ^r)Po - ^-if^a.I + ^l-^- 1 P 0 "■-^i 4 Po + (^a _ ^^a.U 

( 1 iPo 2 " I M Po+ M a) ( l 3 Po 2 -I N Po-N p ) +M q N r p Q 2 +M (1 Z a N 3 


( 9 )- 


Although <Xg S and P s6 are perhaps the most, important of the steady- 
state quantities, three other useful relationships can be derived from 
equations (l), (2), and (5). These are , respectively: 


r 


ss 


Pq 0 ^ S P S S + ^pPo + ^ 

1-Y r 


^o^ss**^a a S6"^ 

I 5^s r ss“ I e p o^s“^^ss" I ^ss" Il p p o 

^ Q 

SS qSb 


C Z r 6r 


( 10 ) 


Character of Results 


As noted, the equations for the quantities of interest are algebraic 
equations which are functions of aerodynamic and inertia parameters and 
the rolling velocity p Q . With known values of the parameters it is 
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possible to assume various values of p Q and compute the corresponding 
magnitudes of <x SB , 3 SS > etc. An Iteration procedure must be used if 

C n is not zero or if the other aerodynamic derivatives are functions 
°a 

of a or p, An illustrative example of the type of information which is 
obtained from such computations is given in figure 3 . In figure 3(a) are 
shown three areas of rolling velocity and the basic factors which govern 
the magnitudes of cxgg and J3 SS in these regions. Of the three, the 
intermediate range of rolling velocity will be of most interest, since 
It is here that the largest changes in Og S and P ss occur. As previously 
mentioned, it is postulated that these changes in magnitude from trim con- 
ditions determine the degree of inertia coupling encountered by an aircraft 
doing aileron rolls with fixed rudder and stabilizer. 

Also shown in figure 3(b) are the two distinctly different types of 
p Q versus 8 a ss plots which can be obtained and which depend primarily on 
whether favorable or unfavorable yaw occurs . These plots are of impor- 
tance since they indicate the magnitude of steady rolling velocity the 
aircraft can achieve. Seme interesting points on roll stability will be 
covered in a subsequent portion of this paper. 


EVALUATION OF THE STEADY-STATE CONCEPT 


To evaluate the method, the steady-state solutions were compared with 
the results of analog- computer studies and with flight data. Although 
these comparisons showed that it was not possible to predict the peak 
magnitudes of the transient a and '(3 with sufficient accuracy for loads 
analysis, the steady-state solution did correctly Indicate the onset and 
severity of roll coupling. A few selected examples follow to illustrate 
the degree of comparison and the usefulness of the steady-state concept. 

The characteristics of the F-100A and the F-1Q2A airplanes were used 
in the calculations and a summary of their mass, geometric, and aerodynamic 
properties is given in tables I and H. However, -unless otherwise noted, 
the basic characteristics of the large- tail F- 10 GA airplane at M = 0.7 

and 30 >000 feet with Cq =0 were used in the following comparisons and 

8 a 

will be identified in the text as the "example aircraft." The steady-state 
a, and p were computed from the complete equations ( 6 ) and ( 7 ) unless 
otherwise noted. 


Comparison of the Steady-State With Analog- Computer Results 


Throughout this section of the report, five-degree- of -freedom (V =0) 
analog-computer results are presented for approximate 3 ^ 0 ° aileron rolls. 
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These 36O 0 rolls are made with pilot applied fixed elevator (or stabilizer) 
and no rudder input during the roll. In figure 4 are shown typical records 
for such a roll with the maximum excursions and average roll rate indi- 
cated. On subsequent figures the maximum excursions in sideslip and in 
angle of attack plus the trim a are used for the comparison with the 
results of the steady- state method. 

Relationship of pitching and yawing frequencies . - References 3, 8, 
and others have shown that it is desirable to have the pitching and yawing 
frequencies of the aircraft equal or close to each other. This conclusion 
is drawn from a study of the equation defining the divergent boundaries , 
a Q = 0, as follows: 

Consider the denominator Sq of the simplified expression for ag B 
or 3 SS (eq. (8) or ( 9 )), 


,+ o + " . + 2 + + , ( ^ 2' f *' \ 

(I1P0 -r^Po+Ma) (isPo -%Po-V+ nq N rPo + Mq z a N p 


( 11 ) 


Now when the normal signs shown above the parameters are examined, it is 
seen that a 0 will always remain positive (and therefore stable) if the 
bracketted terms reverse sign at the same rolling velocity. The rela- 
tionship of the parameters which are necessary for this reversal is seen 
to be: 


_ 3-m 


> 


% 

-Np 




•J-Xe^e 


> (12) 


lY 

Ma 


I Y" I X 


-Np 


T Z 

I Z 

J 



The only manner In which the above relationship can be realized is for 
Iy = Ig and Mq, = -Nj 3 , the latter being the equality between the square 
of the pitching natural frequency and the square of the yawing natural 
frequency. Since 1% is always larger than Iy, the foregoing rela- 
tionship can be approached but not exactly achieved. 
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For later use, the rates of roll at which either of the bracketted 
terms in equation (ll) become zero are defined as critical rolling veloc- 
ities. They are; 


P, 


crit. 


a 


21 n 


J X e w e 

2 ^ I Z -I X^ 


/ -Mq 

J Iy 


( 13 ) 


P crit 


P 


I N ± 


n/i^+^NpIs 

2*3 


I z e UJe 

2( I y - I x ) 



(14) 


Analysis has indicated that when the vertical tail size (i .e . , 



of an airplane is proportioned so as to approximate the best relationship 
between and Np (eq. (12)), excursions in a, and p are alleviated con 

siderably. Since M these results also agree with flight experience, it is 
therefore of interest to examine the effect of arbitrary variations in 
natural frequency parameters M a and Np on the steady-state excursions 
in a and p, and to compare them with analog- computer results . Figure 5 
is presented for discussing the trends which may be expected. Three 
curves have been calculated for both <x B8 and p ss using constants and 
derivatives of the example aircraft, except as noted, starting from 
trimmed level flight. Ctne curve, shown by the solid line, is for the 
best value of Np defined by equation (12) corresponding to the existing 
value of Mq, for the airplane. For this case, Np = -0.75 Mq, and the 
corresponding value of Cn^ is about 15 percent larger than the value 
tabulated in table X for the large- tail airplane. When Np is reduced 

below the best value relative to Hx, a and p excursions are increased. 
In fact, if Np is sufficiently small, instability occurs, as can be 
seen from the dashed curve in figure 5 • The dashed curve, identified as 
Np being equal to 3/10 of the best value relative to Mix, can be asso- 
ciated with the small-tailed F-10QA airplane at the same flight condition 
The remaining curve does not represent a particular configuration of the 
F-100A airplane, but is included to show the effect of increasing Np 
and by a factor of 3 , at the same time retaining the best relation- 
ship between these quantities. Also indicated on the figure for the 
solid and dashed curves are the critical roll rates obtained from equa- 
tions ( 13 ) and (l4) . The critical roll rates for the dash-dot curve are 
out of the range shown. 
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The data of figure 5 clearly show that large roll coupling can occur 
when the best relationship formula between Ng and Mo, is not observed as 
well as the beneficial effects of increasing both the natural frequencies 

in yaw and pitch. ... . -- ... 

In figure 6, the three curves of figure 5 are compared with results 
from five-degree-of-freedam analog- computer solutions. The comparison is 
made in figure 6(a) for the case of Kg having the best value with respect 
to Mo,; in figure 6(b) for Np = 0.3 of the best value; and in figure 6(c) 
for a threefold increase in both Np and while the best-value relation- 
ship is retained. In each of figures 6(a) , 6(b), and 6(c), three curves 
are shown. The. solid line is the calculated steady state, the dashed line 
is the maximum transient excursion for a 36 O 0 roll occurring during the 
departure from trim, and the dash-dot line is the maximum transient excur- 
sion during the return to trim. With the single exception of (3 SS magni- 
tudes at high roll rates (-167° to - 195 ° per second) for the case of 
Np= -0.3 Ma(ls/li) in figure 6(b) , the magnitudes of the maximum a and 3 

transients during the departure frcm trim follow the trends predicted by 
the steady-state equations; those during the return to trim, of course, 
have no steady-state counterpart. It is apparent, however, thht at roll 
rates less then the lower of the two critical roll rates these transients 
appear to deviate from the steady-state curve by roughly the same amount 
as the departure transients, but In an opposite sense. Above this criti- 
cal. roll rate, the return transients follow no obvious pattern with respect 
to the steady-state values . Thus the trends of the maximum angle of attack 
and sideslip as affected by roll rate are correctly given by the steady- 
state method, but the prediction of the magnitudes of these quantities with 
sufficient accuracy for loads analysis is not possible. 

In figure 6(b) (and some subsequent figures) is noted a region in 
which the airplane continued to roll even though the aileron deflection 
was reduced to zero. This phenomenon is associated with favorable yaw 
and possibly could be avoided if the aileron were completely reversed 
rather than returned to neutral. 

Effects of initial angle of attack and elevator deflection .- Previous 
studies (ref. 8) have shown that the initial angle of attack at the onset 
of rolling and inadvertent elevator inputs during a roll can have a large 
influence on the magnitude of the a and 3 excursions during a rolling 
maneuver. It is known further that the excursions will be near a minimum 
if the initial angle of attack is such that the principal, axis and flight 
path are aligned. These effects can be demonstrated also by the use of 
the steady-state method in. the following manner. 

In figure 7 are plotted the ag 6 and 3 SS for Initial angles of attack 
of -4°, +1°, and +6° for the example aircraft which correspond to initial 
angles between the principal axis and the flight path of -5°, 0°, and +5°, 
respectively. As can be noted from the figure, large deviations in ag 8 
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and p ss occur in the region of - 80 ° to -l80° per second roll rate for the 
two cases where the principal axis and flight path are not aligned. Also 
of interest is the reversal of Og S for roll rates greater than -120° per 
second. Since the amount of elevator deflection is directly associated 
■with the initial angle of attack , this change would appear to a pilot roll- 
ing in this region as a control reversal. The large difference in the 
extremes of <Xg S also indicates that the elevator control will be sensi- 
tive in this region. Such characteristics have been observed in flight. 

Figures 8(a), 8(b), and 8(c) have been prepared to show the comparison 
between the steady-state values and the analog- computer results. As is 
seen from these data, the steady-state method is useful in predicting the 
effects of changes in initial angle of attack. 

As an interesting sidelight, the fact that the minimum change in 
Ogg occurs when the principal axis and flight path are nearly aligned 
can be shown easily from the simplified steady-state equation for Og S 
as follows: 

If the damping and engine inertia terms are neglected, the equation 
for Ogg becomes 

ft- J XZ a A 
I g Pp 2 -M M V* % P ° / 

“* B _ (±_ 2 \ 

V-M,, I Y P o J 


When p Q = 0 in this equation, ag S = a trim = orcier then 
to maintain csg S = invariant with roll velocity the terms within 
the parentheses must be equal, or 


1 

M I Y 




Iv 



-1 


This equality reduces to 





T-Xz / declination of the principal 
e , the above equation mean a that if the angle 


Since by definition 
axis below the body axis. 
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of attack of the principal axis is zero, no deviation in ctes with p Q 
■will occur. This conclusion, of course, is modified slightly if all the ... 
terms affecting a^ s are considered. 

Augmentation devices . - The use— of augmenters to reduce the roll- 
coupling tendencies of an aircraft has received considerable attention 
recently (refs. 9 and- 10). The purpose here is to demonstrate that the 
steady-state method provides a means of estimating the effectiveness of 
an augmenter and the magnitude of the control required. Two devices 
have been selected for study; one, the artificial damper in yaw or pitch, 
and the other, a pq plus a pr feedback which tends to cancel certain 
important inertia- coupling terms in the yawing and pitching equations of 
motion, respectively. A simplified block diagram of their mechanization 
for study of five -degree- of -freedom simulation is shown in figure 9 - 

Artificial dampers: Some analog-computer studies (ref. 8 ) have 

indicated that an increase in pitch damping is more effective than an 
equivalent increase in yaw damping in reducing inertia coupling in a 
rolling maneuver. A scrutiny of the simplified steady-state equations 
for ass and p 6S shows that this conclusion should always be true in 
the absence of yawing moments applied by the aileron or rudder (N = 0) . 

For example, increases in pitch and yaw damping by artificial means may 
be approximated in the steady-state equations by increasing the parameters 
Mq and N r , respectively. Referring to equations ( 8 ) and (97 and letting 
N be zero, one sees that Mq occurs twi ce and only id the denominator , 
whereas N r occurs but once in both denominators and in the numerator of 
3 se . Since increases in Mq or N r increase the denominator or numerators 
if the aircraft has not exceeded the divergent boundary, it is apparent 
that increased Mq (pitch damping) will he more effective than N r (yaw 
damping) in reducing either p ss or a s8 . This anomaly is confirmed by 
the data shown in figure 10 , on which have been plotted ags and- Pss ^ or 
the normal example aircraft and for separate increases in the damping in 
pitch and yaw. The increases in the damping ratios noted were computed 
from two- degree- of -freedom analyses in pitch and in yaw at zero roll rate 
and the changes in magnitude were made by appropriate increases in Mq 
and N r (AMq, AW^) in the steady- state equations. Also shown in the figure 
are the steady-state rudder and stabilizer deflections necessary for the 
augmentation as obtained from the following equations: 


■ ss 


AN* 

nT 


* ss 


Ai4 


■'SS 


Ml 4 - 


^6 


The large difference in. the deflections of the control surfaces is due 
primarily to the differences in the control effectiveness. To illustrate 
that the steady-state method will indicate the effect of artificial dampers 
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on roll coupling , figures 11(a ) } 11(b ) f and 11(c) have "been prepared. 

The data in these figures show a good correspondence , trendwise, between 
the steady-state and analog- computer results. For the analog simulation 
the feedback loops C and D of figure 9 were used, with the gains G 1 and 
G 2 adjusted to give the required damping ratios. 

The pq,pr augmenter : It has been shown that an aircraft experiences 

small roll- coupling effects if it rolls at rates well below the critical 
roll velocities defined in equations (13) and (l4). One means then of 
reducing such effects is to increase the critical roll rates by the use 
of an augmenter. As can be seen from the equations defining these criti- 
cal rates, they may be increased in the limit to infinity by canceling, 
with pr and pq feedback, the inertia- coupling terms -I^r and I 3 pq 
in the pitching and yawing equations of motion, respectively (see refs. -9 
and 10 for previous studies on this type augmenter) . Figure 12 illustrates 
the effect of this cancellation on certain steady-state quantities. The 
steady-state rudder and horizontal stabilizer deflections necessary for the 
augmentation were computed from the formulas 



As can be seen from the data in the figure, the pq,pr feedback device 
results in a significant decrease in ag S and 3 SS , but the 5 r ss neces- 
sary to achieve these improvements is large. In an effort to reduce the 
5 r and still retain the desired small variation of ccgg and p ss with 

p Q , artificial damping was added with the favorable results shown in 
figure 13 . 

Figure 9, loops A and B, illustrates the method of simulating the 
pq,pr device on the analog computer. In a practical case the products 
of pr and pq can be obtained by driving the primaries of two transducers 
on the roll rate gyro with amplified signals from the yaw and pitch rate 
gyros . For the analog simulation the gains G x and G 2 were adjusted to 
give exact cancellation of the inertia- coupling terms when the dynamics 
of the servo were neglected. Figure l4 shows the comparison of the steady- 
state and the f ive-degree-of-freedom analog- computer solutions for the two 
types of augmenters . In general, the data in figure 14- along with that of 
figure 11 (a) show that the steady-state method can be used in a qualitative 
manner to assess the roll-coupling effects of these special types of 
stability augmenters. 
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Ability of tie ailerons to create rolling velocity ,- Since roll rate 
has been recognized as a fundamental parameter in the inertia- coupling 
phenomenon, the foregoing evaluation of the steady- state method has been 
predicated on a comparison of the calculated steady-state a and p as a 
function of roll rate with solutions from the five-degree-of -freedom 
analog simulation. However, it is of importance also in the evaluation 
to determine the degree to which the steady-state approach can estimate 
the ability of the aileron to create rolling velocity. In figure 15 are 
shown such rolling performance data. The corresponding comparisons of 
a and p have been given in figure 11(a) . 

The data shown in figure 15 are interesting for several reasons. 
First, the steady-state solution is seen to have a lower and an upper 
statically stable branch (solid lines) connected by a statically unstable 
branch (dash line) . Second, is the close correlation in magnitude between 
the analog-computer' results and those of the lower steady- state branch up 
to about 17° control deflection and the good agreement in trend with the 
upper stable steady-state branch above 30-° aileron-,, deflection. Plotted 
also in the figure is the straight line of the single degree of freedom 
relationship between p Q an! 8^ . . Xb. difference between the straight- 

line and the steady- state solution is due primarily to rolling moments 
induced by sideslip. 

During the course. of obtaining the analog- computer data, seme runs 
were made for rolls greater than 3^0° . These data matched the steady- 
state lower branch closely up to a 6 a of about 25 ° at which point a 
jump in roll rate occurred so. that for larger control deflections the 
steady-state upper branch was followed. Some understanding of the reason 
for this discontinuity in the roll rate was obtaine&Hay ~ considering small 
perturbations about the steady- state values and solving for the chsfracter- 
istic modes of the five equations of motion. The method of reference 6 
was used but only the most important derivatives and inertia terms were 
retained for these small-perturbation computations (l 2 , I 4 , I s , Yp , Yp, 

Y r , J^, Mp, N^, Np, and Lj, were set equal to zero). The results are 

given in figure 16 for the lower statically "stable _ ste‘ady-'state branch. 
These data show that the yaw or dutch roll mode becomes dynamically 
unstable at 27° control deflection. Whether this good agreement between 
the perturbation analysis and the analog- computer results will occur for 
different aerodynamic derivatives is not Imown, as trjansient effects are 
important in determining the critical aileron deflection. 

Since it has been noted by others that the derivative C n has 

5 a 

some influence on the roll-coupling problem, the steady-state solutions 
were computed to illustrate the effect of a change in this derivative at 
both a subsonic and a supersonic speed. These results are given in fig- 
ure 17 and show that a favorable yawing moment due to control deflection 
tends to add favorable yaw (or reduce unfavorable yaw) for the lower 
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branch at either Mach number and so increase the roll rate for a given 
aileron deflection. This effect is so pronounced at the supersonic 
speed that the steady-state solution becomes statically unstable at a 
roll rate of -I 9 O 0 per second with a moderate amount of favorable yawing 
moment. The stability in this instance was checked also by the pertur- 
bation technique at roll rates of -I 9 O 0 and -210° per second. The results 
of these computations are noted in figure 17 and show that at the lower 
roll rate the roll subsidence root is convergent but at - 210 ° per second 
it is divergent. 

The usefulness of the steady-state method and the perturbation 
technique for determining the ability of the aileron to create rolling 
velocity is that it offers a means of assessing the areas of aileron 
deflection (or roll rate) where serious roll coupling occurs . Thus, it 
offers a way of judging where to limit the control deflection (or roll 
rate) if this means is used to avoid roll- coupling difficulties . 


Comparison of Steady-State Solution With Pilot Opinion 


The main objective of the steady-state approach is to provide a 
simple qualitative method whereby in prfeliminary design one can quickly 
assess the onset and degree of roll-coupling effects for any selected 
flight condition. The method would be of questionable value, however, 
if its results were in conflict with pilot opinion regarding actual air- 
plane performance. Seme unpublished pilot- opinion data have been obtained 
from the NACA High-Speed Plight Station at Edwards Air Force Base, 
California. These data, which are for 360 ° rudder-fixed aileron rolls, 
are compared with steady-state solutions in figure 18. Three categories 
of pilot acceptability are shown: intolerable (fig. 18 (a)) ; good 

(fig. 18 (b)); and marginal (fig. 18 (c)). 

In the "intolerable" plot (fig. 18(a)) are the F-102A airplane at 
one flight condition (curve A) and the F-10QA small- tail airplane at two 
flight conditions (curves B and C) . These "intolerable" experiences are 
h.i 1 characterized by having excessively large calculated steady-state (3 
excursions over a sizable roll-rate region attainable by the airplane. 

The good flight experience (fig. 18(b)) have been with the large- tail 
P-100A airplane at two high-speed flight conditions (M = 1.20 and 1.26) , 
and are shown by curves D and E in the figure. These good experiences 
are characterized by small steady-state excursions in the roll-rate region 
attainable. The short vertical lines intersecting the calculated steady- 
state curves designate the average roll rate experienced during the par- 
ticular flight test. For the two conditions shown, these roll rates are 
well below the critical rolling frequency. The marginal experiences 
(fig. 18(c)) have been with the large-tail F-100A airplane at three sub- 
sonic flight conditions (M = 0.70 and 0.93 and altitudes of 30*000 
and 4-0,000 feet). Here, the excursions are relatively large (curves F 
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and C3-) , but these do not caver a broad range of roll rates in the roll-rate 
region attainable by the airplane. It should be noted that for curve G, 
the roll rate experienced was above the critical. For the flight condition 
at M = 0.93 a-t 5-0,000 feet (curve H) , the maximum roll rate experienced 
during the maneuver was on the border of the critical roll-rate region. 

For this condition, the pilot was able to detect an impending undesirable 
condition and so evaluated the flight experience as marginal. 

In summary, one can conclude that pilot opinion is consistent with 
the steady-state concept in these important respects; all divergent con- 
ditions given by steady-state calculations and within the ro lling capa- 
bilities of the aircraft have been recognized as unacceptable by the 
pilots; all. good conditions have been those experienced at high speeds 
where the critical rolling frequencies are high, and at roll rates which 
have been appreciably below these critical frequencies. All marginally 
acceptable conditions have been those in which roll rate either has 
remained just below the critical regime or has passed unobtrusively 
through it . 

These marginal cases could very well become unacceptable for roll 
maneuvers initiated fr.cm initial conditions other than those shown, since 
the aircraft is rolling at rates in the critical region where initial 
angle of attack, inadvertent elevator motions', and other effects can 
cause an appreciable change in the magnitude of the actual excursion 
from trim. 


CONCHJSIONS 


An analytical method has been developed herein for use in analyzing 
the roll-coupling phenomenon. This method is based on the c oncept that 
variations in the calculated steady-state angle of attack and sideslip 
with roll rate give a measure of the onset and degree of roll coupling. 
From the results of comparisons of the steady-state values with those 
from five-degree-of-freedom analog- computer studies and from limited 
flight data the following conclusions are drawn: 

1. The trends of the maximum angles of attack and sideslip as 
affected by roll rate are correctly given by the steady-state method, 
but the prediction of the magnitudes of these quantities with sufficient 
accuracy for loads analysis is not possible. 

2. An analysis of the steady-state equations provides an under- 
standing of and a useful guide for determining the effects of various 
aerodynamic and inertia parameters on roll coupling. 

3 . The method appears useful to define critical flight conditions 
and to assess the merits of various stability augmenters devised to 
minimize roll- coupling effect; 
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The method, because of the ease and rapidity with which 
calculations can be made, appears particularly promising for preliminary 
design work in defining areas in which a complete simulator study should 
be undertaken . 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif. , Nov. 7 , 1 956 
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APPENDIX A 

THE COMPLETE . COEFFICIENTS OP THE. STEADY- STATE 
SOLUTION FOR Ogg AND p 8S 

The constants a Q , t! 0 , and c Q of equations (6) and (7) are: 

a o = I xisPo 4 i)Po 3+ l (l _ Y r ) Po 2+ 

[ (1-Y r ) (l 3 MpZ a -I^)+I M N p -Z a Yp(l 1 I N +l3l M )+MpN r ]p 0+ 

{[ ( 1- Y r ) N p +N r Y p ] ( -l^+Mq^) +I N Z a [ I M Y p - (l-Y r )M p ]j- 

b Q = [I 2 I 3 (l-Y r )-I 1 I 3 Yp]p 0 4 +[-I 2 I K (l-Y r )-I 1 l 3 Y+Yp(l 1 I F +l 3 l M ) ]p 0 3 + 

[ (1-Y r ) (-I 3 M-I 2 Np-M a N p )+Y(l 1 I N +l 3 l M )+Y p (l 1 N p -I M I N -M < 1 N r ) + 

Yp (-I 2 N r -IiN p -I 2 .I 3 Z) ]p 0 a +[ (1-Y r ) (jff lf -NM q _-IgZMp-N p Mp)+Y(l 1 N 3 -I M I N -M ( 1 N r )+ 
Y p ( -IwNp-NrMp ) +Y p ( -I a.N+1 ^Z+IsImZ+I^p ) ] p Q + [ (l-Y r ) (MNp-M^ZNp+IflMpZ- 
NMp ) +Y ( - Ij^p-N^ ) +Y p (MN r -M q ZN r -I M I w Z+I M N) ] 

c Q = |[I 2 I 3 (l-Y r ) -IiI 3 Yp ] Z^+I i(l 3 Z+N p ) +I 2 N r j- P 0 3 + 1[ -I 2 I N ( 1- Y r ) -I 1 I 3 Y+ 

Ypd^M+IaI M ) ]Z a - I M( I 3 z+N p) +I i( N - I N Z )}p o 2+ t ( 1_Y r ) (-Is^+M^p+Is^Z- 
MgZaFp ) +ZqY ( I iIjj+I 3 I^) +Yp( -I^I^Z a -M q NrZ a ,+I^xNr) “ (MNr - MqZNp-Ij ( jI^Z+ 

I m N) ]p Q +[ (l“Y r ) (l^-M q Z a N+I^ a -I I ^Z)+Y(-I M I N Z a -M q N r Z a +H x N r ) ] 
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TABLE I.- SUMMARY QF PHYSICAL PROPERTIES AND DERIVATIVES 
FOR THE F-lOQA AIRPLANE 
[All derivatives not listed are taken as zero] 

Physical Properties 


c 

11.33 


11,050 

J xz 

508 

b 

36.58 

% 

59,000 

H 

X 

n> 

17,550 

S 

376.0 

I Z 

67,000 

m 

7^5 


Derivatives 



M = 0.64 

M = 0.70 

M = 0.93 

M = 1.20 

M = 1.26 


-3.89 

-4.27 

-4.87 

-3.84 

-3.84 

c Dit 

— 

-1.00 . 

— 

- .60 

— 


1 -3.60 

-3.75 

-5-55 

-4.15 

-4.15 


-.29 

-.42 

-.81 

-1.22 

-.91 

Sx 

— 

-I.25 

— 

-1.35 

— 

Cn P 

0 

-.025 

0 

- - 

0 

0 

Cn r 

-.16 

-.30 

-.29 

-.19 

-.19 

(large tail) 

— 

.095 

.109 

.097 

.086 

Cnp (small tail) 

.039 

.04-3 

— 

— 

— 

°X 


-.006 

— 

-.0095 

— 

Cn & 

°r 

— 

-.0315 

— 

— 

— 

Cy 

• Y P 

0 

.17 

0 

0 

0 

C Y r 

0 

• 34 

0 

0 

0 

% 

0 

-.62 

0 

0 

0 

C lp 

— 

-.29 

— 

-.406 

— 

C lr 

— 

.12 

— 

• 03 

1 — 

c lp at a = 0 

— 

-.008 

— 

-.050 

— 


— 

-.516 

— 

-.332 

— 

V 

i 

-.044 

— 

-.038 

— 
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TABIE II.- SUMMARY OF PHYSICAL PROPERTIES AND DERIVATIVES 
FOR THE F-102A AIRPLANE AT M = O.T5 
[All derivatives not listed, are taken as zero] 


Physical Properties 


c 

23.1 


13,600 

I xz 

2310 

b 

37.8 

Iy 

89 >00 

x X e w e 

17,550 

S 

661.5 


99,600 

m 

711 


Derivatives 


Ax 

-2.81 


-1.20 


-.258 

Cn r 

-.069 

On 

^3 

.065 


r* 
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Primary influence on magnitude 




Figure 3 .- Character of results of steady-state computations. 
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Figure 4.- Time histories of typical 3^0° roll maneuver showing values 
of a and 0 used for comparing five- degree- aT -freedom analog- 
computer results with steady- state calculations. 
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0 -40 -80 -120 -160 -200 

p 0 , deg/sec 


Figure J Effect of nonalignment of principal and flight-path axes for 
the example aircraft using the simplified steady-state equations. 
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Steady state 

Maximum excursions during 

— departure from trim 

— return to trim 


Pav or Po> deg /sec 


(a) a -€ = 0 

trim . .. 

Figure 8.- Comparison of steady-state results with five-degree- of 
freedom analog- computer results for the example aircraft with 
different alignments of principal and flight-path axes. 
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Maximum excursions during: 



0 -40 -80 -120 -160 -200 

P 0 V or Po . deg/sec 

(c) a trim -e * ' 5 -°° 


Figure 8.- Concluded. 
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Figure 9 .- Simplified block diagram of a stability augmenter incorporating damping and artificial 
nullification of inertia cross-coupling terms (pq.,pr feedback) . 
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Maximum excursions during : 


departure from trim 



O -40 -80 -120 -160 -200 

P av or Po> deg/sec 

(a) Without augmentation; elevator and rudder 
fixed; £ a = 0 . 25 ; = 0.11. 

Figure 11.- Comparison of steady-state results with five-degree-of -freedom 
analog- computer results for the example aircraft with different amounts 
of damping in pitch and yaw. 
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°0 -40 -80 120 160 200 

P QV or p o> deg/sec 

(b) Pitch, damping increased to £ a = 0.5; £p = 0.11. 
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Figure 11.- Continued. 
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Figure 12.- Comparison of i 
aircraft with one augmei 
inertia- coupling terms. 
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(a) The pq,pr cLevi 

Figure l4.- Comparison of steady-state resull 
analog -computer results for the example a" 
Mlity augmenters. 
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(b) The pq,pr device plus added damping; t = f = 0 

W K 

Figure l4.- Concluded. 
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-100 


-200 


F-IQ2 A 

A M = 0.75 ; h p = 39,500 ; I g 
F- 100 A, smoll tail 

B M = 0.70 ; h p = 32,000 ; Ig 
C M ■ 0.64 ; h p = 10,000 ; 2.8 g 




200 0 


p. , deg/sec 


Not divergent^ 
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-200 


F - 100 A, large toil ; I g F- IOOA, large tail ; I g 

D M B 1.26 ; h p = 40,000 ft F M = 070 j h p = 40,000 ft 

E M = 1.20 ; h p = 30,000 ft 6 M s 0.70; h p * 30,000 ft 

H M = 093; h p = 40,000 ft 


(a) Intolerable. 


(b) Good. 


(c) Marginal. 


Figure 18.- CanrpariBon of steady- state solution with pilot opinion for several airplanes; 

360° rudder- fixed aileron rolls. 
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